The fluctuations, under time reversal, of the natural time and the entropy distinguish 
similar looking electric signals of different dynamics 
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We show that the scale dependence of the fluctuations of the natural time itself under time reversal 
provides a useful tool for the discrimination of seismic electric signals (critical dynamics) from noises 
emitted from manmade sources as well as for the determination of the scaling exponent. We present 
recent data of electric signals detected at the Earth's surface, which confirm that the value of the 
entropy in natural time as well as its value under time reversal are smaller than that of the entropy 
of a "uniform" distribution. 
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I. INTRODUCTION 

In a time series comprising TV events, the natural time 
Xk = k/N serves as an index^ 3 - for the occurrence of 
the fc-th event. In natural time analysis, the time evolu- 
tion of the pair of the two quantities (xfc, Qk) is con- 
sidered, where Qk denotes in general a quantity pro- 
portional to the energy released during the fc-th event. 
In the case of dichotomous electric signals (e.g., seis- 
mic electric signal (SES) activities, i.e., low frequency 
< 1Hz electric signals that precede earthquakes, e.g., see 
RefsiiiMlilllilllil Q k stands for the duration of 
the fc-th pulse (cf. The SES activities should not be con- 
fused with pulses of very short durations observed some 
minutes before earthquakes^) . It has been showni£ that 
natural time domain is optimal for enhancing the signals' 
localization in the time-frequency space, thus conforming 
to the desire to reduce uncertainty and extract signal in- 
formation as much as possible. The entropy S in natural 
time is defined^ as the derivative with respect to q of 
the fluctuation function (x q ) — (x) q at g = 1: 



S= (xlnx) - (x) ln(x> 



(1) 



where (/(%)) = Y,k=iPkf(Xk) and p k = Qk/J2n=iQn- 
It is dynamic entrop y 17 ' 18 and exhibits^ concavity, pos- 
itivity and Lesch o 20 ' 21 stability. Note that S should not 
be confused with Cov[x, lnx] = (x m x) — (x)( m x) since 
in general (lux) m (x)- The value of the entropy upon 
considering the time reversal T, i.e., Tpk = Pn— fe+ii is 
labelled by SL . The value of 5_ i a 19 ' 22 i 23 , in general, dif- 
ferent from S, and thus S does satisfy the conditions to be 
"causal" in the following sense (see RefGJ and references 
therein) : When studying a dynamical system evolving in 
time, the "causality" of an operator describing this evo- 
lution assures that the values assumed by the operator, 
at each time instant, depends solely on the past values of 
the system. Hence, a "causal" operator should be able 
to represent the evolution of the system according to the 
(true) time arrow, thus the operator can represent a real 
physical system evolving in time and reveal the differ- 
ences arising upon time-reversal of the series. 

The statistical properties of S and S- have been stud- 



ied in a variety of model s 22 ' 23 . In the case of a "uniform" 
distribution S = SL = S u (= In 2/2 - 1/4 w 0.0966). 
The "uniform" distribution (defined in RefsfHl2~4l) has 
been analytically studied in and corresponds to 

the case when Qk are independent and identically dis- 
tributed (IID) positive random variables of finite variance 
including the case of Markovian dichotomous electric sig- 
nals studied in Reffl6l. The "uniform" distribution cor- 
responds to p(x) = lj where p{x) is a continuous proba- 
bility density function (PDF) corresponding to the point 
probabilities pk used so far. When Qk of a "uniform" dis- 
tribution are perturbed by a small linear trend, we find 22 
(see also Eq.@, below) that (S - S U )(S- - S u ) < (cf. 
this simple example, which shows that S captures the ef- 
fect of a linear trend, may be considered as clarifying the 
meaning of S, see Section V of Refl22T). Another model 
studied is when the increments of Qk are positive IID, 
in this case we find 23 that S « 0.048 and SL w 0.088 
which are both smaller than S u . The same holds, i.e., 
that both S and S— are smaller than S u , in the examples 
of an on-off intermittency model discussed in Reft22l as 
well as for a multiplicative cascades modelS 3 . adjusted to 
describe turbulence data. 

A case of practical importance is that of the SES activ- 
ities. SES activities (critical dynamics) exhibit infinitely 
ranged long-range temporal correlation a 22 ' 23 ' 24 which are 
destroyed 2 ^ after shuffling the durations Qk randomly. 
An interesting property emerged from the data analysis 
of several SES activities refers to the fact^ 2 . that both S 
and S- values are smaller than the value of S Ul i.e., 



(2) 



in addition to the fact that for SES activities^! 3 - the 
variance 



— /,,2 



(X 2 ) - (X) 2 » 0.070. 



(3) 



These findings -which do not holdi* for "artificial" 
noises (AN) (i.e., electric signals emitted from manmade 
sources)- have been supported by numerical simulations 
in fractional Brownian motion (fBm) time serie o 22 ' that 
have an exponent cxdfa, resulted from the Detrended 
Fluctuation Analysis (DFA) 2 ^^, close to unity. This 
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model have been applied since fBm (with a self-similarity 
index H m 1) has been found^I as an appropriate type 
of modeling process for the SES activities. These sim- 
ulations resulted in values of S and S— that do obey 
relation (2) (see Fig.4 of ReflU) and m ~ 0.070 (see 
Fig. 3 of Refl23l). It was then conjectured^ that the va- 
lidity of the relation (2) stems from infinitely ranged 
long-range temporal correlations (cf. H ~ 1). On the 
other hand, for short-range temporal correlations (e.g. 
when modeling Qk by an autoregressive process Qk = 
aQk-i + 5fc + c, \a\ < 1 and c stands for an appropriate 
constant to ensure positivity of Qk or Qk = \aQk-i + 9k\ 
where gk is Gaussian IID variables) the values of both 
S and S- approach (see Appendix A) that of S u and 
Kl — > k«j where k u = 1/12 denotes the corresponding 
value of the "uniform" distribution^. 

The scope of this paper is twofold: First, in Section 
II, we point out the usefulness of the study of the fluc- 
tuations of the natural time itself under time reversal. 
In particular, it enables the determination of the scaling 
exponent, thus allowing the distinction of SES activities 
from similar looking AN. Second, in Section III, we pro- 
vide the most recent experimental data that strengthen 
the validity of the relations (2) and (3) for SES activities. 
The earthquakes that followed the latter SES activities 
are described in Section IV. Section V, summarizes our 
conclusions. 



II. THE FLUCTUATIONS OF NATURAL TIME 
UNDER TIME REVERSAL 

The way through which the entropy in natural time 
captures the influence of the effect of a small linear trend 



has been studied, as mentioned, in RefH2 on the basis 
of the parametric family of PDFs: p(x', e) = 1 + e(x — 
1/2), where e measures the extent of the linear trend. 
Such a family of PDFs shares the interesting property 
7~p{x] e) = p(x;~ e )i i- e ; the action of the time rever- 
sal is obtained by simply changing the sign of e. It 
has been shown^ that the entropy 5(e) = S\p(x',e)], as 
well as that of the entropy under time reversal S- (e) = 
S[Tp(x;e)], S-(e) — S(—e), depend non-linearly on the 
trend parameter e: 




However, it would be extremely useful to obtain a linear 
measure of e in natural time. Actually, this is simply the 
average of the natural time itself: 



f 1 e 

(x) = J o xp(x; t)dx = 2 + ^2 ■ ( 5 ) 



If we consider the fluctuations of this simple measure 
upon time-reversal, we can obtain information on the 
long-range dependence of Qk- We shall show that a 
measure of the long-range dependence emerges in nat- 
ural time if we study the dependence of its fluctuations 
under time-reversal Axf = E[((x) — (Tx)) 2 } on the win- 
dow length I that is used for the calculation. Since 
Tp k =Pi-k+i, we have 



A x f = E[« X ) - (T X )) 2 ] = E 



T (Pk-Pl-k+l) 



(6) 



where the symbol E[. . .] denotes the expectation value obtained when a window of length I is sliding through the 
time series Qk- E[. . .] is well defined when all the {pk}k=i,2,...i involved in its argument are also well defined. The 
evaluation of E[. . .] can be carried out either by full or by Monte Carlo calculation. In order to achieve this goal, from 
the original time-series {Qk}k=i,2,...L, we select segments {<9m -M-i}i=i,2,..J of length i, and the argument of E[. . .] 
is computed by substituting pk — Q ma +k-i/ 5Zi=i Qmo+i-l- The sum of the resulting values over the number of the 
selected segments (different mo) is assigned to E[. . .]. The full calculation refers to the case when mo takes all the 
L — I + l(mo = 1, 2, . . . L — / + 1) possible values, whereas the Monte Carlo when mo is selected randomly. 
By expanding the square in the last part of Eq. (|6|) , we obtain 



A Xl = X! ( T ) ^^P k -Pl-k+l) 2 } + X! ~JT E ^P k -Pl-k+l)(Pm -Pl-m+1 



k=l 



k^m 



(7) 



The basic relation^ that interrelates pk is Pk = 1 or equivalently pk = 1 — Jm/fcP" 1, ^ subtracting from the 
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last expression its value for k = I — k + 1, we obtain pk — pi-k+i = — J2 m ^k(P™ ~ Pi-m+i), and thus 

(Pk - Pi-k+i) 2 = —/ l(j>k - Pi-k+i)(p m -pi- m +i). (8) 

By substituting Eq.® into Eq.©, we obtain 

^ f \ ^ ^ J^jYYIj 

A Xl = ~2J ( T ) /_/ E t( Pfc - Pl-k+l)(Pm ~ Pl-m+l)} + 22 ~JF E [(P fe -Pl-k+l)(Pm ~ Pl-m+l)} (9) 



which simplifies to 

Ax? = — -jp^- E[(p fc -pz-fc+i)(p m -^-m+i)] 

fe,m 

(10) 

The negative sign appears because (pfc — pi-k+i) and 
(Pm — Pi-m+i) are in general anti-correlated due to 
Eq.||5J). Equation ([TO]) implies that A^ 2 measures 
the long-range correlations in Qk'. If we assume that 

-E[(p k - Pl-k+l)(p,n ~ Pl-m+l)} OC (k — m) 2xH /P (cf. 

Pk scales as l/l, e.g. seei£), we have that 

A x f oc l 4+2 x»/l 4 (11) 

so that 

A X i(=^/a^)^I xh , (12) 
where is a scaling exponent. 



A. Fractional Brownian motion and fractional 
Gaussian noise time series 

In order to examine the validity of the above result 
Eq. (|12[) when Qk are coming from fBm or fractional 
Gaussian noise (fGn), we employed the following pro- 
cedure: First, we generated fBm (or fGn) time-series 
Xk (consisting of 2 x 10 4 points) for a given value of 
H using the Mandelbrot- Weierstrass functio n 28 i 29 i 30 as 
described in Refl22l. Second, since Qk should be posi- 
tive, we normalized the resulting time-series to zero 
mean and unit standard deviation and then added to the 
normalized time-series Nk a constant factor c to ensure 
the positivity of Qk = Nk + c (for the purpose of the 
present study we used c = 10). The resulting Qk time- 
series were then analyzed and the fluctuations of Axi 
versus the scale I are shown in Figs. QJa) and QJd) for 
fGn and fBm, respectively. The upper three panels of 
FigO] correspond to fGn while the lower three to fBm. 
We observe (see Fig[TJb)) that for fGn we have the in- 
terconnection: xh ~ H — 1 corresponding to descending 
curves(see FigJIJa)), whereas for fBm the interconnec- 
tion turns (see FigOJe)) t° : Xh ~ H corresponding to 
ascending curves(see Figfljd)). 



In order to judge the merits or demerits of the proce- 
dure proposed here for the determination of the scaling 
exponent, we compare FigsOJb) and[lje) with FigsOJc) 
andQJf), respectively, that have been obtained by the 
well-established DFA method 25 - 26 . This comparison re- 
veals that the results are more or less comparable for fGn, 
while for fBm the exponent xh deviates less from the be- 
havior of an ideal estimator of the true scaling exponent 
(drawn in dashed green) compared to cldfa> especially 
for the largest H values. 



B. The fluctuations of the natural time to 
distinguish seismic electric signal activities from 
similar looking AN 

The physical meaning of the present analysis was fur- 
ther investigated by performing the same procedure in 
the time-series of the durations of those signals analyzed 
in Refll9i that have enough number of pulses e.g.« 10 2 
(cf. the signals depicted in FigJ3]could not be analyzed in 
view of the small number of pulses) . The relevant results 
are shown in Fig[5J Their inspection interestingly indi- 
cates that all seven AN correspond to descending Axi 
curves versus the scale I, while the three SES activities 
to ascending curves (in a similar fashion as in FigsQJa) 
and [IJd) , respectively) as expected from the fact that 
the latter exhibit^ 4 - infinitely ranged long-ranged tempo- 
ral correlations (having H close to unity), while the for- 
mer do not. Hence, the method proposed here enables 
the detection of long-range correlations even for datasets 
of small size (« 10 2 ), thus allowing the distinction of SES 
activities from AN. 



III. RECENT DATA OF SEISMIC ELECTRIC 
SIGNALS ACTIVITIES 

First, Fig|3Ja) depicts an electric signal, consisting of 
a number of pulses, that has been recorded on Novem- 
ber 14, 2006 at a station labelled^ PIR lying in western 
Greece (close to Pirgos city). This signal has been clearly 
collected at eleven measuring electric dipoles with elec- 
trodes installed at sites that are depicted in a map given 
in Ref 31. The signal is presented (continuous line in 
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red) in Fig. 02 a ) in normalized units, i.e., by subtract- 
ing the mean value and dividing by the standard devi- 
ation. For the reader's convenience, the corresponding 
dichotomous representation is also drawn in Fig. 02a) 
with a dotted (blue) line, while in Fig. 03c) we show (in 
red crosses) how the signal is read in natural time. The 
computation of S and S— leads to the following values: 
S = 0.070 ± 0.012, S- = 0.051 ± 0.010. As for the vari- 
ance «i, the resulting value is «i = 0.062 ± 0.010. These 
values more or less obey the conditions (2) and (3) that 
have been found to hold for other SES activities 22 . Note 
that the feature of this SES activity, it is similar to the 
one observed at the same station before the magnitude 
M « 6.7 earthquake that occurred on Jan 8, 2006, see 
ReflH. 

A closer inspection of Fig. 02 a ) reveals the following 
experimental fact: An additional electric signal has been 
also detected (in the gray shaded area of Fig. 02 a ))j 
which consists of pulses with markedly smaller ampli- 
tude than those of the SES activity discussed in the pre- 
vious paragraph. This is reproduced (continuous line in 
red) in Fig. 02^) in an expanded time scale and for the 
sake of the reader's convenience its dichotomous repre- 
sentation is also marked by the dotted (blue) line, which 
leads to the natural time representation shown (dotted 
blue) in Fig. 02 c )- The computation of S and S- gives 
S = 0.077 ± 0.004, S- = 0.082 ± 0.004, while Kl is found 
to be ki = 0.076 ± 0.005. Hence, these values also obey 
the conditions (2) and (3) for the classification of this 
signal as an SES activity. 

The two aforementioned signals have been followed by 
two significant earthquakes as described in Section IV. 
This conforms to their classification as SES activities, 
which has been completed in an early version of this 
paper— on November 16, 2006. 

Second, very recently, i.e., on July 2, 2007 and July 
10,2007, two separate electric signals were recorded at a 
station labelled PAT lying in central Greece (close to Pa- 
tras city) at 38.32°N21.90°E. The signals are presented 
(continuous line in red) in Fig. 02(d) andOH(e) in normal- 
ized units in a similar fashion as in Figs02 a )j(b). Their 
corresponding dichotomous representation are also drawn 
with dotted (blue) lines, while in Fig02f) we show (in 
red crosses and blue asterisks, respectively) how the sig- 
nals are read in natural time. The computation of Ki, 
S and S- leads to the following values: For the signal 
on July 2, 2007: ki = 0.072 ± 0.005, S = 0.073 ± 0.007, 
S- = 0.081 ± 0.006, for the signal on July 10, 2007: 
Ki = 0.073±0.004, S = 0.085±0.005, SL = 0.080±0.004. 
An inspection of these values reveals that they obey the 
conditions (2) and (3) and hence both signals can be clas- 
sified as SES activities. The procedure for the current 
study of the subsequent seismicity that occurred after 
these SES activities is described in the next Section. 



IV. THE SEISMIC ACTIVITY THAT 
FOLLOWED THE SES ACTIVITIES 

We discriminate that during the last decade SES activ- 
ities are publicized only when their amplitude indicates 
that the impending earthquake has an expecte d 13 ' 23 mag- 
nitude comparable to 6.0 unit or larger. 



A. The case of the SES activities of Figs. 3(a), (b) 

According to the Athens observatory (the data of 
which will be used here), a strong earthquake (EQ) with 
magnitude 5.8-units occurred at 13:43 UT on February 3, 
2007, with epicenter at 35.8°N 22.6°E, i.e., almost 80 km 
to the southwest of the 6.9 EQ of January 8, 2006, (cf. the 
magnitude announced from Athens observatory is equal 
to ML+0.5, where ML stands for the local magnitude). 
This was preceded by a 5.2-units EQ that occurred at 
22:25 UT on January 18, 2007 at 34.8°N 22.7°E. The 
occurrence of these two EQs confirm the classification 
as SES activities of the signals depicted in Figs02 a ) an d 
02b). (Note that preseismic information based on SES ac- 
tivities is issued only when the magnitude of the strongest 
EQ of the impending EQ activity is estimated -by means 
of the SES amplitude 3 ^- to be comparable to 6.0 units or 
larger—.) 

Here, we show that the occurrence times of the afore- 
mentioned two EQs can be estimated by following the 
procedure described in RefsfH l22ll23ll35l and u sing the or- 
der parameter of seismicity proposed in Ref[35|, i-e, the 
normalized power spectrum in natural time n(<^) as 
<f> ^ (see also below). We study how the seismic- 
ity evolved after the recording of the SES activities on 
November 14, 2006, at PIR station (which were classi- 
fied as SES activities in the initially submitted version 
of the present paper on November 16, 2006). The study 
is made either in the area A^j^'lf E^gf or in the area 
^NljUgElf;^ (see FigH]), by considering three magni- 
tude thresholds M t hres = 3.2, 3.4 and 3.6 (hence six 
combinations were studied in total). If we set the natu- 
ral time for seismicity zero at the initiation of the SES 
activity at 17:19 UT on November 14, 2006, we form 
time series of seismic events 3 -^ in natural time for various 
time windows as the number N of consecutive (small) 
EQs increases. We then compute the normalized power 



spectrum 



1.22,23.35 



in natural time H(4>) for each of the 



time windows. Excerpts of these results, which refer to 
the values during the periods: (a) December 25, 2006, 
to January 17, 2007, and (b): January 18 to January 
31, 2007 are depicted with red crosses in FigfSJ respec- 
tively. This figure corresponds to the small area B with 
Mthres — 3.4. In the same figure, we plot in blue the 
normalized power spectrum obeying the relation^' 2 -' 3 -' 3 ^ 



HH 



18 
5^ 



6cosa-> 12 sin w 



5lu 3 



(13) 



which holds when the system enters the critical 
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stage (u> — 2tt(/), where (f> stands for the natural 
frequencji 2 -! 3 -^). The date and the time of the occur- 
rence of each small earthquake (with M > 3.4) that oc- 
curred in the area B , is written in green in each panel (see 
also Table IJ). An inspection of FigfS^a) reveals that the 
red line approaches the blue line as N increases and a co- 
incidence occurs at the small event of magnitude 3.7 that 
occurred at 03:22 UT on January 17, 2007, i.e., roughly 
two days before the 5.2-units EQ at 22:25 UT on Jan- 
uary 18,2007. A similar behavior is observed in FigfSJb) 
in which we see that a coincidence occurs at the small 
event of magnitude 3.6 at 18:40 UT on January 31, 2007, 
i.e., roughly three days before the strong EQ of magni- 
tude 5.8-units that occurred at 13:43 UT on February 3, 
2007. To ensure that these two coincidences in FigsfSJa) 
and (b) are true ones^ ' 13 i 31 ' 35 (see also below) we also 
calculate the evolution of the quantities ki,S and S- and 
the results are depicted in FigJSjfor the three magnitude 
thresholds for each of the aforementioned two areas A 
and B. 

The conditions for a coincidence to be considered as 
true are the following (e.g., see RefQ], see also 3 ' 13 ' 31 ' 35 ): 
First, the 'average' distance (D) between the empirical 
and the theoretical 11(0) (i.e., the red and the blue line, 
respectively, in FigfSJ) should be^ ' 13 ' 31 ' 35 smaller than 
or equal to 10~ 2 . See FigJTJ where we plot (D) versus 
the conventional time during the whole period after the 
recording of the SES activities on November 14, 2006, for 
both areas, i.e., the large one (area A) and the small (area 
B) and the three magnitude thresholds. For the sake of 
the readers convenience, the mean value of the results 
obtained for the three thresholds is also shown in black. 
Second, in the examples observed to date^i 3 - ' 13 ' 31 ' 34 ' 35 , a 
few events before the coincidence leading to the strong 
EQ, the evolving 11(0) has been found to approach that 
of Eq.(l), i.e., the blue one in Figf5] , from below (cf. 
this reflects that during this approach the Ki-value de- 
creases as the number of events increases). In addition, 
both values S and SL should be smaller than S u at the 
coincidence. Finally, since the process concerned is self- 
similar (critical dynamics), the time of the occurrence 
of the (true) coincidence should not change, in princi- 
ple, upon changing either the (surrounding) area or the 
magnitude threshold used in the calculation. Note that 
in FigJTl at the last small events ,i.e., the rightmost in 
Figs|U[a) and 03b), respectively (i.e., the magnitude 3.7 
event on January 17, 2007 and the second event of magni- 
tude 3.6 on January 31, 2007) just before the occurrences 
of the 5.2-units and 5.8-units EQs, in both areas A and 
B, the mean value (see the black thick lines in FigJT]) of 
(D) obtained from the three magnitude thresholds be- 
come smaller than or equal to 10 -2 . Hence, these two 
coincidences can be considered as true. 



In summary, the SES activities recorded on November 
14, 2006, at PIR station (presented in Figs(21^a),(b)) have 
been followed by two EQs with magnitudes 5.2-units and 
5.8-units that occurred on January 18 and February 3, 
2007. The time of the occurrences of these two EQs are 
determined within a narrow range of a few days upon 
analyzing, in natural time, the seismicity subsequent to 
the SES activities. 



B. The case of the SES activities of Figs. 3(d), (e) 

The actual amplitude (in mV) of the most recent SES 
activities recorded at PAT on July 2, 2007 and July 10, 
2007 (see Fig|3]Jd) and (e), respectively) can be visual- 
ized in Figs(5fa) and[5fb) where the original recordings 
of a measuring electric dipole (with length L «5km) are 
reproduced. For the sake of comparison, in Figs(5{c),(d) 
we also present the corresponding SES activities at the 
same station, i.e., PAT, that precede d 23 ' 31 the magni- 
tude 6.0-class earthquakes that occurred with epicenters 
at 37.6°N20.9°E on April 11 and 12,2006. Furthermore, 
in FiglSJe), we show the SES activity 3 ! at PAT on Febru- 
ary 8, 2007, which was followed by a magnitude class 6.0 
earthquake at 38.3°N20.4°E that occurred on March 25, 
2007. 

In order to determine the occurrence time of the im- 
pending EQs, we currently apply the procedure explained 
in the previous subsection by studying the seismicity in 
the areas A, B, C (see Fig|H|). Since the result should 
exhibit spatial scale invariance, the epicenter (s) will lie 
either in the area B or in C depending on whether the 
areas A and B or A and C show true coincidence. 
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FIG. 1: (color online) Examples of log-log plot of the fluctuations Ax; of the natural time under time reversal versus the scale I 
for fGn (a) and fBm (d). (b) and (e) depict the values of the scaling exponent \h, introduced here, versus the true exponent H 
for fGn and fBm, respectively. For the sake of comparison, (c) and (f) are similar to (b) and (e), respectively, but for the DFA 
exponent an fa- The (blue) dotted curves show the ±cr deviation from the average value (obtained after 10 2 runs) depicted 
by the (red) solid curves. The (green) dashed straight lines correspond to the ideal behavior of each exponent and have been 
drawn as a guide to the eye. 
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FIG. 2: (color online) The log-log plot of A\i versus the scale I for three SES activities (Kl, K2 and U) and seven AN (nl-n5, 
n7 and nlO) treated in RefQjJ (cf. these signals have enough number of pulses in order to apply the present analysis). The 
values of A\i are divided by the corresponding values A\5 at the scale / = 5. 



V. CONCLUSIONS 

First, the scale dependence of the fluctuations of the 
natural time under time reversal distinguish similar look- 
ing electric signals emitted from systems of different dy- 
namics providing a useful tool for the determination of 
the scaling exponent. In particular, SES activities (criti- 
cal dynamics) are distinguished from noises emitted from 
man-made electrical sources. 

Second, recent data of SES activities are presented 



which confirm that the value of the entropy in natural 
time as well as its value under time reversal are smaller 
than that of a "uniform" distribution. 
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FIG. 3: (color online) (a) The electric signal recorded on November 14, 2006 at PIR station (sampling rate f exp — 1Hz). The 
signal is presented here in normalized units (see the text). The corresponding dichotomous representation is shown with the 
dotted (blue) line. The gray shaded area shows an additional signal (consisting of pulses of smaller duration) superimposed 
on the previously mentioned signal, (b) The signal belonging to the gray shaded area in (a) is given here in an expanded time 
scale, while its dichotomous representation is marked by the dotted (blue) line, (c) How the signals in (a) and (b) are read 
in natural time: the continuous (red crosses) and dotted (blue asterisks) bars correspond to the durations of the dichotomous 
representations marked in (a) and (b), respectively. 



APPENDIX A: THE CASE OF SIGNALS THAT 
EXHIBIT SHORT-RANGE TEMPORAL 
CORRELATIONS 

Here, we present results of modeling Qk by short- 
ranged temporal correlated time-series. Two examples 
were treated by numerical simulation: (i) A stationary 
autoregressive process Qk — aQk-i + 9k + c, \a\ < 1, 



where gk are Gaussian IID random variables, and c 
stands for an appropriate constant to ensure positivity 
of Q k . (ii) Qk = \aQk-i + 9k\- FigureEJa) depicts the 
results for S,S— and K\ for the first example versus the 
number of Qfc, whereas FigfTUTb) refers to the second ex- 
ample. In both cases S and S— converge to S u whereas K\ 
to the value k u — 1/12 corresponding to the "uniform" 
distribution. For the reader's convenience, the values of 
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FIG. 4: Map of the area surrounding the station PIR (solid dot) and the epicenters of the EQs (stars) with magnitude 5.2-units 
and 5.8-units that occurred on January 18 and February 3, 2007, respectively. The seismicity subsequent to the SES activities 
on November 14, 2006, has been studied in the gray shaded areas (the large and the small area are designated A and B, 
respectively) . 




FIG. 5: (color online) The normalized power spectrum(red) H((f>) of the seismicity as it evolves event by event (whose date 
and time of occurrence are written in green in each panel) after the initiation of the SES activity on November 14, 2006. The 
two excerpts presented here, refer to the periods: (a) December 25, 2006, to January 17, 2007, and (b): January 18 to January 
31, 2007, and correspond to the area B with M t hres = 3.4. In each panel only the spectrum in the range 4> € [0,0.5] (for the 
reasons discussed in Refs l 35-) is depicted (separated by the vertical lines), whereas the 11(0) of Eq, (113|l is depicted by blue 
color. The minor horizontal ticks for (f> are marked every 0.1. 
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FIG. 6: (color online) Evolution of the quantities k\, S and 5*_; they are shown in (a), (b) and (c), respectively, for three 
magnitude thresholds, i.e., M > 3.2, 3.4 and M > 3.6, for both areas, i.e., the large (area A) and the small (area B). The solid 
horizontal lines correspond to ki = 0.07 and S u = ln2/2 — 1/4. 



S u and k u are designated by the horizontal solid black APPENDIX B: THE SEISMIC ACTIVITY THAT 
lines. FOLLOWED THE RECENT SES ACTIVITIES OF 

FIGS.8(A),(B) 



Considering the Athens observatory preliminary cata- 
logue, the seismic activity (M t hres > 3.0) that occurred 
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FIG. 7: (color online)The average distance (D) between the calculated and the theoretical H(<t>) curves versus the conventional 
time. The calculation of (D) is made upon the occurrence of every consecutive earthquake when starting the calculation after 
the initiation of the SES activities on November 14, 2006 (depicted in Figs|3la),(b)) for the large (area A) and the small (area 
B) for the three magnitude thresholds Mthres = 3.2, 3.4 and 3.6. The black solid line corresponds to the mean value obtained 
from the three thresholds. 



in area A (see Fig. 9) after the initiation of the SES activ- 
ity on July 2, 2007 (Fig.8(a)) until 03:27 UT of September 
25, 2007 is shown in Fig fTTT a). The evolution of the cor- 
responding parameters K\, S, S- and (D) calculated for 
three magnitude thresholds, i.e.,M t ^ res = 3.0,3.1 and 3.2 
are shown in Figs lllf b). (c) and (d) respectively. To in- 
vestigate the spatial invariance, the computation was re- 
peated for several smaller areas, three of which are shown 
in Figs (HUB! and [14] (which are different from the areas 
B and C of Fig. 9) along with the evolution of the corre- 
sponding parameters. The same was repeated for an area 
(see Fig fTS)) somewhat larger than A. An inspection of all 
these figures, i.e., FigsfTT1tdT5l suggests that presumably 
a true coincidence has just been approached, thus being 
very close to the critical point. 



APPENDIX C: THE SEISMIC ACTIVITY THAT 
FOLLOWED THE MOST RECENT SES 
ACTIVITIES AT PAT AND PIR 

Here, we report the update results of the seismic ac- 
tivity that followed the SES activity at PAT on October 
9, 200822 and the SES activity at PIR on December 12, 
2008^2 by following the procedure described by Sarlis et 
al4£. The subsequent seismicity of the former SES activ- 
ity was studied in the area Ngy fE^g g while that of the 
latter in the selectivity map of PIR depicted in Fig. [TH] 
The results, when considering the seismicity until early in 
the morning of February 2, 2009, for magnitude thresh- 
old M t /i res =3.3, are shown in Figs. [T7] and [TS] for the 
former and the latter SES activities at PAT and PIR, re- 
spectively. An inspection of these figures reveals that in 
both areas the probability Prob(Ki) versus K\ -calculated 
in all the possible regions of each area as described in 
Ref4£- maximizes at n\ w 0.070 upon the occurrence of 
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FIG. 8: Recent SES activities recorded at PAT. 



the events marked with arrows, thus probably indicating 
the approach to the critical point. 

Note added on February 20, 2009: Actually, at 23:16 
UT on February 16, 2009 a strong earthquake with 
magnitude Ms(ATH)=6.0 (i.e., M L (ATH)=5.5) occurred 
with an epicenter at 37.1°N20.8°E, which clearly lies in- 
side the selectivity map of PIR depicted in FigfTBl 

We now present the results, when considering the seis- 
micity until early in the morning of February 19, 2009. 
Figure [T91 depicts the results for Mt^ res =3.2 in the se- 



lectivity map of PAT (i.e., in the area N^fE^;!), which 
show that the aforementioned probability Prob(Ki) ver- 
sus K± maximizes at k± ~ 0.070 upon the occurrence of 
a M L =3.2 event at 00:35 UT on February 19, 2009 with 
epicenter at 37.8°N21.2°E. As for the results in the selec- 
tivity map of PIR (i.e., the shaded area in Fig fTr?)) . they 
are shown in Figj20] for M t / lres =3.4 and reveal the fol- 
lowing: After the aforementioned strong earthquake on 
February 16, 2009, the bimodal curve Prob(rei) versus k\ 
shows a clear local maximum at n\ « 0.070 upon the oc- 
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FIG. 9: A map showing the areas A, B and C at which the analysis of the seismicity in natural time is currently carried out. 
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FIG. 10: (color online) The thick colored lines in (a) and (b) depict the average value of S (red),.?- (blue) and ki (green) 
versus the number of Qk for the two examples (i) and (ii), respectively, mentioned in the Appendix. The thinner colored lines 
refer to the ±cr deviation from the average value. 



currence of the Ml=3.5 event at 00:30 UT on February 
18, 2009 with an epicenter at 34.9°N23.4°E. This local 
maximum does exhibit magnitude threshold invariance 
(since it also appears when investigating M t ^ res =3.3 and 
M thres =3.5). 

Updated note on March 13, 2009: Several hours be- 
fore the occurrence of the aforementioned magnitude 6.0 
earthquake on February 16, 2009, an electrical anomaly 
of significant amplitude, see Fig Ell appeared at PIR. It 
basically consists of three bay like transient changes, thus 
having a form essentially different than the SES activi- 
ties that preceded the major earthquakes on February 14, 
2008 and June 8, 2008. This electrical anomaly could be 



in principle attributed to the 6.0 earthquake on Febru- 
ary 16, 2009 that occurred almost 11 hours later, but 
this earthquake was also preceded by the SES activity at 
PIR on December 12, 2008. Alternatively, this anomaly 
might be related to a new impending earthquake. For 
this reason, a natural time analysis of the seismicity sub- 
sequent to the electrical anomaly of February 16, 2009 
was carried out. In this calculation we considered the 
events that occurred in the PIR selectivity map shown 
in Fig[22] (cf. this, which is slightly different than that 
depicted in Fig. 16, has been drawn by taking into ac- 
count the whole region to the west of the Hellenides). 
The results of the analysis shown in Fig|2"3l indicate that 



13 



for M t / lres =3.6 and 3.5 (see a and b, respectively) the 
probability Prob(Ki) maximizes at k x « 0.070 upon the 
occurrence of a Ml=3.6 event at 02:48 UT on March 9, 
2009, while for M t / lres =3.3 the maximization occurs upon 
the occurrence of a Ml=3.3 event on March 10, 2009. 
Interestingly, a similar analysis (subsequent to the SES 
activity at PAT) in the selectivity map of PAT, i.e., in the 



area Np-gE^g - !, reveals that a maximization of Prob(Ki) 
at Kj. as 0.070 also occurs on March 10, 2009 (see FiglMj). 
We are currently investigating the spatial and magnitude 
threshold invariance of the aforementioned maxima to ex- 
amine whether they actually point to an approach to the 
critical point. 
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FIG. 11: (color online) Results ol the analysis of seismicity in the dark gray shaded area of (a) after the SES initiation on July 2, 2007. 
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FIG. 12: (color online) Results of the analysis of seismicity in the dark gray shaded area of (a) after the SES initiation on July 2, 2007. 
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FIG. 13: (color online) Results of the analysis of seismicity in the dark gray shaded area of (a) after the SES initiation on July 2, 2007. 
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FIG. 14: (color online) Results of the analysis of seismicity in the dark gray shaded area of (a) after the SES initiation on July 2, 2007. 
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FIG. 15: (color online) Results of the analysis of seismicity in the dark gray shaded area of (a) after the SES initiation on July 2, 2007. 
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FIG. 16: The selectivity map of PIR (shaded area) along with 
the selectivity map of PAT (rectangular area) 
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78 3.7 2009/2/2 



0:42:44 
8:54:47 
23:40:43 
18:58:05 
0:57:07 
1:28:55 



N(deg)E(deg) 
38.14 20.35 
37.79 21.27 
37.78 21.09 
38.05 20.42 
38.26 20.42 
38.51 20.42 
37.59 20.89 
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FIG. 17: The Prob(Ki) versus Ki of the seismicity for 
M t h res =3.3 subsequent to the SES activity recorded at PAT 
on October 9, 2008 within the selectivity map of PAT shown 
in Fig. 1161 For the sake of clarity, only the last 7 events are 
depicted. The arrows mark the maxima at ki =0.070 (see the 
text) 
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No MLDate(Y/M/D)Time(UT) N(deg)E(deg) 
87 3.3 2009/1/27 0:50:31 36.59 21.70 
8 3.3 2009/1/27 19:32:10 36.60 21.75 

89 3.3 2009/1/27 23:40:43 38.05 20.42 

90 3.3 2009/1/28 18:58:05 38.26 20.42 

91 3.3 2009/1/29 01:15:04 36.59 21.79 

92 3.4 2009/1/30 14:49:36 37.01 20.90 

93 3.3 2009/1/31 15:38:54 36.55 21.67 

94 3.4 2009/1/31 20:24:26 35.70 22.10 

95 3.3 2009/1/31 20:29:03 35.77 23.15 

96 4.8 2009/2/1 0:57:07 38.51 20.42 

97 3.3 2009/2/1 6:16:31 36.14 21.50 

98 3.7 2009/2/2 1:28:55 37.59 20, 




FIG. 18: The Prob(«i) versus K\ of the seismicity for 
M t h res =3.3 subsequent to the SES activity recorded at PIR on 
December 12, 2008 within the selectivity map of PIR shown 
in Fig. 1161 For the sake of clarity, only the last 12 events are 
depicted. The arrows mark the maxima at ki =0.070 (see the 
text). 



No ML Date(Y/M/D)Time(UT)N(deg)E(deg) 

125 3.7 2009/2/2 01:28:55 37.59 20.89 

126 3.9 2009/2/2 04:13:34 37.61 20.95 

127 3.2 2009/2/4 08:38:41 38.20 21.84 

128 3.2 2009/2/14 08:49:23 37.79 21.14 

129 3.2 2009/2/18 08:51:16 37.56 22.77 

1 30 3.2 2009/2/1 9 00:35:55 37.83 21 .22 




FIG. 19: Prob(Ki) versus k\ of the seismicity for Mth res =3.2 
subsequent to the SES activity recorded at PAT on October 9, 
2008 within the area N||;fE?9;|. For the sake of clarity, only 
the last 6 events are depicted. Prob(fti) exhibits a maximum 
at ki=0.070 upon the occurrence of the 3.2 event on February 
19, 2009 (see the text). 
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FIG. 20: Prob(fi:i) versus m of the seismicity for M t h res =3.4 
subsequent to the SES activity recorded at PIR on December 
12, 2008 within the selectivity map of PIR shown in Fig. 1161 
For the sake of clarity, only the last 17 events are depicted. 
The arrows mark the maxima of Prob(Ki) at ki =0.070 (see 
the text). 
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FIG. 21: Electrical recordings from PIR; for the sites of the 
electrodes see Ref.— . 




20°E 22°E 24°E 26°E 28°E 



FIG. 22: The selectivity map of PIR when including the whole 
region to the west of Hellenides (shaded area) . The rectangu- 
lar area corresponds to the selectivity map of PAT. 




FIG. 23: Prob(Ki) versus Ki of the seismicity in the shaded 
region of Fie(22lfor M*k...>«=3.6 (a). 3.5 (b) and 3.3 (c) subse- 




FIG. 24: Prob(Ki) versus K\ of the seismicity for Mt hres =3.3 
subsequent to the SES activity recorded at PAT on October 9, 
2008 within the area N^fE?!;! . For the sake of clarity, only 
the last 6 events are depicted. Prob(Ki) exhibits a maximum 
at K\ =0.070 upon the occurrence of the event at 15:22 UT on 
March 10, 2009. 
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TABLE I: The catalogue (Mthres = 3.2) of the Institute of Geodynamics of the National Observatory of Athens (GI-NOA) 
(available from http://www.gein.noa.gr/services/monthly-list.html on February 5, 2007) for the Large area (area A) 
during the following period: From the initiation of the SES activity almost at 17:00 UT on November 14, 2006 (Fig.l of the 
main text) until the occurrence of the 5.8-units earthquake on February 3, 2007. Note that Ms(ATH)=M+0.5, where M stands 
for the local magnitude ML or the "duration" magnitude MD. Two events reported by GI-NOA: at 17:32:23.7 on January 28, 
2007, and at 5:11:06.5 on January 30, 2007, were omitted since they could not be independently confirmed by other seismic 
networks operating at that time. 



N 


Year 


Month 


Day 


UT 


Lat.(°N) 


Lon.(°E) 


depth(km) 


M 


1 


2006 


11 


16 


4:20:40.0 


35.66 


23.26 


54 


3.5 


2 


2006 


11 


16 


12:30:50.3 


35.44 


23.42 


5 


3.7 


3 


2006 


11 


16 


18:50:02.5 


36.97 


23.11 


25 


3.7 


4 


2006 


11 


16 


21:53:12.5 


35.72 


24.33 


10 


3.4 


5 


2006 


11 


17 


16:31:26.4 


36.69 


23.32 


24 


3.2 


6 


2006 


11 


18 


3:48:58.3 


34.79 


24.07 


16 


3.6 


7 


2006 


11 


18 


21:51:21.8 


34.78 


24.36 


34 


3.2 


8 


2006 


11 


19 


9:31:09.1 


35.73 


21.95 


19 


3.6 


9 


2006 


11 


19 


10:54:52.3 


35.99 


23.65 


10 


3.2 


10 


2006 


11 


19 


14:11:45.4 


35.75 


22.55 


24 


3.3 


11 


2006 


11 


21 


5:57:05.5 


34.81 


24.28 


54 


3.3 


12 


2006 


11 


21 


8:42:06.0 


37.42 


21.96 


19 


3.2 


13 


2006 


11 


21 


12:50:34.0 


35.45 


23.90 


10 


3.3 


14 


2006 


11 


24 


17:16:58.4 


37.16 


21.58 


12 


4.1 


15 


2006 


11 


24 


17:42:56.3 


35.60 


23.50 


7 


3.3 


16 


2006 


11 


24 


20:26:02.9 


37.12 


21.56 


11 


3.5 


17 


2006 


11 


24 


21:41:44.1 


37.21 


21.63 


8 


3.3 


18 


2006 


11 


25 


3:14:31.2 


34.94 


23.91 


5 


3.2 


19 


2006 


11 


25 


4:55:01.3 


37.15 


21.86 


9 


3.3 


20 


2006 


11 


25 


9:30:18.1 


36.75 


21.73 


12 


3.6 


21 


2006 


11 


26 


17:58:05.5 


35.45 


23.30 


11 


4.1 


22 


2006 


11 


26 


19:45:54.7 


35.40 


23.27 


2 


3.3 


23 


2006 


11 


26 


23:57:49.1 


35.47 


23.36 


5 


3.4 


24 


2006 


11 


27 


3:12:49.8 


35.42 


23.31 


9 


3.3 


25 


2006 


11 


27 


11:58:03.7 


35.53 


23.43 


24 


3.2 


26 


2006 


11 


28 


9:07:38.6 


36.07 


22.35 


10 


3.3 


27 


2006 


11 


28 


15:07:43.1 


36.06 


22.42 


43 


3.6 


28 


2006 


11 


28 


15:42:59.4 


35.71 


22.11 


10 


3.4 


29 


2006 


11 


28 


22:54:24.0 


37.06 


21.51 


12 


3.2 


30 


2006 


11 


29 


14:11:02.8 


35.97 


23.11 


19 


3.5 


31 


2006 


11 


30 


23:30:58.6 


35.43 


23.32 


3 


3.2 


32 


2006 


12 


1 


2:37:02.1 


35.58 


23.60 


10 


3.5 
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TABLE I: Continued 



N 


Year 


Month 


Day 


UT 


Lat.(°N) 


Lon.(°E) 


depth(km) 


M 


33 


2006 


12 


1 


19:57:55.4 


37.07 


21.47 


3 


3.2 


34 


2006 


12 


1 


22:49:35.6 


35.72 


22.52 


114 


3.7 


35 


2006 


12 


2 


0:29:20.8 


34.74 


22.87 


45 


3.2 


36 


2006 


12 


2 


10:30:35.3 


37.42 


22.04 


4 


3.3 


37 


2006 


12 


3 


1:53:22.4 


36.72 


21.81 


24 


3.5 


38 


2006 


12 


3 


2:37:36.7 


36.40 


21.66 


10 


3.2 


39 


2006 


12 


4 


19:32:24.0 


34.99 


23.44 


14 


3.4 


40 


2006 


12 


7 


22:22:49.4 


35.73 


23.17 


32 


3.2 


41 


2006 


12 


8 


10:54:40.0 


34.95 


23.42 


5 


3.3 


42 


2006 


12 


8 


12:48:15.5 


37.02 


21.93 


10 


3.2 


43 


2006 


12 


8 


19:57:27.2 


35.06 


23.50 


7 


3.6 


44 


2006 


12 


9 


22:45:56.5 


35.89 


23.38 


24 


3.4 


45 


2006 


12 


10 


3:29:46.1 


34.84 


24.48 


10 


3.2 


46 


2006 


12 


10 


12:41:27.6 


35.80 


23.04 


59 


3.2 


47 


2006 


12 


11 


3:36:22.0 


35.74 


24.08 


16 


3.4 


48 


2006 


12 


11 


19:11:29.4 


35.93 


23.40 


15 


3.2 


49 


2006 


12 


12 


3:40:55.0 


35.54 


22.81 


17 


3.7 


50 


2006 


12 


13 


9:45:40.0 


37.25 


21.94 


20 


3.2 


51 


2006 


12 


14 


16:44:02.0 


35.95 


23.46 


10 


3.3 


52 


2006 


12 


15 


4:40:12.5 


35.01 


23.37 


80 


3.4 


53 


2006 


12 


15 


10:23:32.1 


36.09 


22.25 


45 


3.7 


54 


2006 


12 


16 


10:43:15.4 


34.85 


24.33 


6 


3.6 


55 


2006 


12 


16 


15:23:33.2 


34.92 


23.45 


5 


3.7 


56 


2006 


12 


17 


4:44:38.7 


34.73 


24.00 


30 


3.5 


57 


2006 


12 


17 


4:44:42.7 


35.15 


24.12 


59 


3.5 


58 


2006 


12 


17 


10:20:49.7 


36.06 


21.76 


37 


3.3 


59 


2006 


12 


17 


14:13:23.4 


36.21 


21.70 


10 


3.4 


60 


2006 


12 


17 


20: 0:20.6 


34.84 


24.21 


34 


4.0 


61 


2006 


12 


17 


22:38:37.5 


36.65 


21.00 


10 


3.2 


62 


2006 


12 


20 


4:34:09.2 


37.48 


21.50 


10 


3.2 


63 


2006 


12 


20 


4:45:15.0 


35.42 


21.43 


10 


3.5 


64 


2006 


12 


20 


19:27:32.5 


34.60 


23.79 


11 


3.8 


65 


2006 


12 


20 


20:17:02.7 


36.68 


21.50 


4 


3.2 


66 


2006 


12 


24 


6:58:02.3 


34.94 


24.03 


10 


3.6 


67 


2006 


12 


24 


7:12:11.4 


36.29 


22.20 


24 


3.8 


68 


2006 


12 


25 


14:08:59.4 


34.83 


22.66 


39 


4.4 


69 


2006 


12 


25 


14:15:50.1 


34.99 


23.04 


39 


4.0 


70 


2006 


12 


25 


14:18:50.7 


35.09 


23.06 


27 


4.0 


71 




1 9 


9^ 




OO.Ul 


99 Qz, 


1 8 




72 


2006 


12 


25 


14:57:00.5 


35.04 


22.84 


14 


4.1 


73 


2006 


12 


25 


17:29:17.9 


35.92 


23.62 


5 


3.2 


74 


2006 


12 


25 


18:13:44.5 


35.09 


23.29 


10 


3.7 


75 


2006 


12 


28 


6:32:19.4 


37.53 


21.81 


10 


3.2 


76 


2006 


12 


29 


2:00:54.9 


36.77 


21.81 


5 


3.3 


77 


2006 


12 


30 


22:51:38.0 


35.39 


23.34 


13 


3.2 



29 



TABLE I: Continued 



N 


Year 


Month 


Day 


UT 


Lat.(°N) 


Lon.(°E) 


depth(km) 


M 


78 


2006 


12 


31 


7:00:51.9 


35.85 


22.08 


10 


3.2 


79 


2006 


12 


31 


18: 6:27.5 


35.14 


22.76 


12 


3.7 


80 


2007 


1 


1 


15:27:48.7 


34.72 


24.12 


23 


3.5 


81 


2007 


1 


3 


15:04:04.9 


35.37 


23.29 


10 


3.4 


82 


2007 


1 


3 


18:18:34.5 


36.54 


21.74 


26 


3.8 


83 


2007 


1 


4 


7:57:08.3 


37.37 


21.56 


15 


3.3 


84 


2007 


1 


4 


14:55:24.6 


36.29 


21.92 


10 


3.2 


85 


2007 


1 


4 


17:42:54.7 


34.91 


23.61 


45 


3.3 


86 


2007 


1 


4 


20:41:40.9 


34.89 


24.13 


48 


3.2 


87 


2007 


1 


4 


22:47:39.1 


36.96 


21.07 


25 


3.2 


88 


2007 


1 


7 


14:07:11.5 


37.10 


21.93 


26 


3.2 


89 


2007 


1 


8 


16:13:02.8 


35.10 


23.03 


10 


3.4 


90 


2007 


1 


9 


9:55:29.2 


36.70 


21.55 


31 


3.4 


91 


2007 


1 


9 


15:54:47.3 


35.89 


23.61 


3 


3.2 


92 


2007 


1 


9 


23:09:20.0 


36.21 


22.70 


116 


3.6 


93 


2007 


1 


11 


5:50:39.5 


35.02 


22.48 


38 


3.8 


94 


2007 


1 


13 


11:12:43.3 


35.48 


23.51 


8 


3.4 


95 


2007 


1 


14 


9:09:23.0 


35.30 


23.38 


3 


3.8 


96 


2007 


1 


14 


16:43:01.6 


35.06 


23.20 


85 


4.1 


97 


2007 


1 


15 


0:55:18.8 


37.47 


21.07 


15 


3.2 


98 


2007 


1 


15 


3:13:45.8 


36.06 


22.40 


13 


3.4 


99 


2007 


1 


15 


6:56:47.0 


35.42 


23.56 


33 


3.3 


100 


2007 


1 


15 


17:50:48.7 


36.56 


21.65 


5 


3.4 


101 


2007 


1 


17 


1:52:18.8 


36.20 


21.58 


34 


3.7 


102 


2007 


1 


17 


3:22:53.9 


35.34 


23.53 


12 


3.7 


103 


2007 


1 


18 


22:25:23.0 


34.84 


22.67 


39 


4.7 


104 


2007 


1 


18 


23:31:32.7 


34.76 


21.36 


10 


3.4 


105 


2007 


1 


19 


8:42:22.8 


34.72 


22.57 


38 


4.4 


106 


2007 


1 


20 


7:05:37.0 


34.95 


24.09 


19 


3.5 


107 


2007 


1 


20 


7:24:24.6 


34.86 


23.96 


19 


3.4 


108 


2007 


1 


20 


10:24:55.4 


36.88 


22.10 


10 


3.2 


109 


2007 


1 


20 


11:57:23.7 


36.69 


22.47 


10 


3.5 


110 


2007 


1 


20 


23:50:04.1 


34.93 


24.40 


28 


3.3 


111 


2007 


1 


21 


17:52:34.9 


35.46 


23.46 


11 


3.4 


112 


2007 


1 


25 


10:24:47.3 


35.33 


23.40 


2 


3.3 


113 


2007 


1 


25 


21:09:45.7 


35.09 


23.25 


22 


4.3 


114 


2007 


1 


28 


15:19:23.6 


36.29 


23.17 


25 


3.4 


115 


2007 


1 


28 


19:12:16.4 


36.95 


21.13 


10 


3.4 


116 


2007 


1 


30 


5:11:08.9 


35.94 


23.41 


28 


3.3 


117 


2007 


1 


30 


10:38:42.7 


36.98 


21.10 


12 


3.5 


118 


2007 


1 


30 


21:11:43.0 


36.13 


22.17 


30 


3.3 


1 1 Q 


9007 


i 
i 




i.ou.oy. ± 


AQ 


99 78 


ou 


3 A 


120 


2007 


1 


31 


14:45:09.7 


34.67 


22.42 


29 


3.6 


121 


2007 


1 


31 


18:40:54.1 


36.26 


22.48 


23 


3.6 


122 


2007 


2 


1 


16:45:35.3 


34.69 


22.46 


21 


3.9 


123 


2007 


2 


1 


21:06:54.7 


36.06 


21.62 


42 


3.2 


124 


2007 


2 


2 


6:39:05.6 


35.21 


23.25 


33 


3.3 


125 


2007 


2 


2 


13:27:53.1 


37.22 


21.63 


4 


3.4 



